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ABSTRACT 

Recent indications from both particle physics and cosmology suggest the existence of more than 
three neutrino species. In cosmological analyses the effects of neutrino mass and number of species 
can in principle be disentangled for fixed cosmological parameters. However, since we do not have 
perfect measurements of the standard A Cold Dark Matter model parameters some correlation re- 
mains between the neutrino mass and number of species, and both parameters should be included in 
the analysis. Combining the newest observations of several cosmological probes (cosmic microwave 
background, large scale structure, expansion rate) we obtain N^s — 3.58^q :[g(68% CL)1;q;53(95% CL) 
and < 0.60eV(95%CL), which are currently the strongest constraints on N^s and from 

an analysis including both parameters. The preference for Ncff> 3 is now at a 2a level. 



1. INTRODUCTION 

Neutrinos are the lightest massive known particles, and 
despite the fact that we know they have mass, they are 
treated as exactly massless by the Standard Model of 
particle physics. Neutrino oscillation experiments us- 
ing solar, atmospheric, and reactor neutrinos have mea- 
sured mass differences between the three Standard Model 
specieij^to be Am^^ = |(2.43lJj;J|)^ IQ-^I eV^ and 



Am^i = (7.50 ± 0.20 ) x 10"^ eV^ ( |Fukuda et al.||l998 
Beringer et al. 2012). The Heidelberg-Moscow experi- 
ment has limited the mass of the electron neutrino to 
be less than 0.35 eV using neutrino-less double /3— decay 
( Klapdor-Kleingrothaus & Krivosheina 2006), but no 
current laboratory experiment has the sensitivity to mea- 
sure the absolute neutrino mass. 

The oscillation data from short baseline experiments 
exhibit some tension allowing for, or even preferring, the 
existence of additional neutrino species. Depending on 
the exact analysis the preferred scenarios include either 
one or two sterile neutri nos in addition to the three nor- 
mal ones (3-f-l or 3-1-2) (Kopp et al. '2011', 
20TT| |Huber|[20Tr ' 



Mention et al. 



Giunti & Laveder^^2011^ 



Ubservations of the Cosmic Microwave Background 
(CMB) also seem to favour additional species of radiation 
to be present at the time of decoupling over and above 
that provided by the photon density and three Standard 
Model neutrinos; most notably the results from the Ata- 
cama Cosmology Telescope (ACT) QDunkley et al.|2011|) 
of iVeff = 5.3 lb 1.3 but oth e r authors e.g. [Komatsu et al.| 
( |20lT| ); |Hou et al.| ([20TT| ; [Keisler et al.| ( I201ip tinTa 
similar preference. N^gis the parametrisation m terms 
of the effective number neutrino species. In analogy to 
dark energy and dark matter, this extra radiation is often 
dubbed dark radiation. 
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^ 'I'echnically they are not Standard Model neutrinos because 
they have mass, but here we refer to the Ue, u^, and Ut as Standard 
Model, and any additional species as sterile. 



The neutrino oscillation results favour a large mass 
difference (e.g. Kopp et al. 2011) providing a lower 
limit on the sterile neutrino mass of the order of 1 eV, 
which is incompatible with cosmological mass constraints 
from combinations of CMB and Large Scale S tructure 
(LSS) measurements (Y] rrii, < 0.3 



0.6 eV, 



S0rensen et al. 2012 Archidiacono et al. 2012 



e.g 



Riemer- 
Joudak: 



|et al.|20T2] ), but it has been shown that factors such as 
initial lepton asymmetry can possibly alleviate these con- 
straints by introducing a non-thermal neutrino velocity 
spectrum and thereby a voiding the LSS mass constraints 
( [Ilannestad et al.||2012[ ). 

Many recent publications have addressed the issues of 
neutrin o mass and effective number of neutrinos sepa- 
rately (iMoresco et al. 2012aj iXia et al.||2012[ |de Putter 



et al. 2O12~Riemer-S0rensen et al. 2012) but very few 



consider the correlation between the two. In this paper 
we address these correlations and demonstrate that due 
to imperfect measurements of the standard A Cold Dark 
Matter (ACDM) parameters, the two parameters are not 
entirely independent but should be addressed simultane- 
ously. Sec. [2 explains thephysical and observable effects 
of iVcff and ^Jtoj/, Sec. [3] and Sec. H describe the data 
and methods used to obtain the results of cosmological 
fits which are presented in Sec. [5] 

2. PHYSICAL AND OBSERVABLE EFFECTS 

This section briefly explains the physical effects of 
J2 n^v and A^eff and how they can be measured. 

The parameter that is probed directly by cosmologi- 
cal analyses is not the neutrino mass, but the neutrino 
density, p^, which can be exp ressed in terms of the mass 
( Lesgourgues k Pastor|[2006] ): 



Pc 



93.14 cV/i2 



(1) 



where pc is the critical energy density for a flat Universe, 
and Ni, is the number of massive neutrino states. Be- 
cause of the smallness of the measured mass differences 
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Figure 1. Illustration of how the CMB and matter power spectra change for varying neutrino mass (solid lines) and effective number of 
neutrinos (dashed lines) fixing all other parameters (to WMAP 7-yr values for ACDM). ^m,^ does not affect the CMB power spectrum 
much, but changes the matter power spectrum significantly. The effect of N^g is clearly visible for small scales (high values of I) in the 
CMB power spectrum, and the two parameters are clearly degenerate in particular for the matter power spectrum. The shaded regions 
indicate the normalised uncertainties of current experiments. 



relative to the upper limits, it is reasonable to assume 
that that if the individual neutrino masses are near the 
upper limit, they are effectively equal. 

2.1. Massive neutrino effects on the expansion 

Due to the weakness of the Weak interaction the neu- 
trinos decouple early from the tightly coupled baryon- 
photon fluid {T = 2 — 4 MeV) , while the expansion is still 
dominated by radiation. If ^ > 1 eV the neutrinos 
become non-relativistic before recombination. However, 
cosmological constraints already indicate ^ < 1 eV 
so the neutrinos become non-relativistic after recombina- 
tion, and while they are relativistic their energy density 
contributes as radiation rather than matter. The pres- 
ence of "neutrino radiation" affects the expansion rate 
and changes the time of the matter-radiation equality. 
The latter determines when the density perturbations 
can begin to collapse gravitationally and structures can 
evolve. The higher ^rui, or A'cff, the higher the ra- 
diation density and the later the equality takes place. 
Changing the expansion rate and time of the matter- 
radiation equality has some measurable consequences de- 
scribed below and illustrated in Fig. [T] 

• The position of the peak of the matter power spec- 
trum is determined by the particle horizon length 
scale at the time of matter-radiation equality. Mov- 
ing the equality to later times shifts the peak to 
larger scales. 

• The expansion rate determines the time available 
for sound waves to propagate in the baryon-photon 
fluid before recombination, which will imprint on 
the BAO scale measurable from both the CMB and 
the large scale structure. Increasing the expansion 
rate will move the BAO peak to smaller scales. 

• Delaying matter-radiation equality enhances the 
ratio of even to odd peaks in the CMB power 



spectrum because the dark matter potentials have 
longer time to grow, and thereby increasing the 
BAO amplitude (the gravitational potentials are 
deeper, but the photon restoring force is the same). 
It also increases the early Sachs- Wolfe effect, which 
enhances the overall normalisation of the CMB 
power spectrum. 

• Changing the expansion rate can significantly af- 
fect the Big Bang Nucleosynthesis (BBN) leading 
to different deuterium and helium abundances than 
the observationally co nfirmed p redictions of stan- 
dard BBN ( .Lcsgourgues fc Pastor 2006; Stcigman] 
20121). 



2.2. Massive neutrino effects on structure formation 

When the neutrinos have de coupled, the ir thermal ve- 
locities decay adiabatically as (Komatsu et al. 20ir|: 

wthormai = 151(1 + ^)(leV/mi,) kmsec"^ . (2) 

When they become non-relativistic they behave as a 
species of warm/hot dark matter, suppressing density 
fluctuations on scales smaller than their free-streaming 
lengt h (at the time when th ey become non-relativistic) 
( jLesgourgues fc Pastor|[2006l ): 



«FS ^ 



Hit) 



iO.82 



2 Wthermal(i)(l + z) 



(3) 



(i^).Mpc-M4) 



Neutrinos cannot cluster on scales smaller than their 
free-streaming length and consequently the density fluc- 
tuations are suppressed by a factor proportional to the 
neutrino density. This is observable in the density power 
spectrum as a suppression of structure on small scales 
(see left part of Fig.[T]). However, the presence of neutri- 
nos also enhances the density perturbations (by adding to 
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the radiation density for fixed flm), which boosts s truc- 
ture forma tion in general, so for higher iVofi, erg (Hou 
et al. 112011 1 will increase 



i'he evolution of the density perturbations can be pre- 
dicted by solving the coupled set of Einstein and Boltz- 
man equations. For small density perturbations the lin- 
earised equations are sufficient (e.g. using CAMBI, but 
when the density perturbations grow sufficiently~large, 
the higher-order, non-linear terms become important and 
have to be included, which is computationally very inten- 
sive (and not yet solved analytically). Instead o ne has 
to rely on second order pert urbation theory (e.g 



Saito 



et al 



et al 



2009[ |Taruya et al.|2012 ) or simulations (e.g. pjmitl 



wm 



Jennings et al.| |2010 
JNumerical solutions show that for = $1^/17^ < 0.07 
the sup pression is dP/P = —8/^ for linear structure for 



mation ( Hu et aL| [l998) and the effect increases for non- 
linear structure formation (8 9.6 Brandbyge et al. 
Brandbygc & Hanncstad 2009 ; Brandby ge et ai" 
Viel et al.||2010(|Agarwal fc Feldman||201lt 



2008 



2010 



The signature of tree-streaming massive neutrinos can 
be detected in both the power spectrum of CMB tem- 
perature anisotropies and in the (matter) density power 
spectrum of large-scale structure. As illustrated in Fig.[T] 
the larger the neutrino mass, the more prominent these 
effects on the power spectra of the CMB and LSS. 

The density power spectrum can be traced either 
through galaxies (galaxy power spectrum) or through 
the absorption lines from neutral hydrogen clouds along 
the line of sight to distant quasars (Lyman-alpha forest 
power spectrum). The neutrinos affect the total mat- 
ter power spectrum, which may differ slightly from the 
tracer power spectrum that we actually measure. De- 
pending on the tracer, the bias, Pobs = b^Pm, might be 



scale dependent or scale-independent ( Schulz & White 
2006| . 

Xnother effect that is important for spectroscopic 
galaxy surveys is the redshift space distortions caused by 
galaxies falling into the gravitational potential of galaxy 
clusters ( |Kaiser|1987{|Percival fc White 2009). As struc- 



ture formation becomes non-linear, so do the redshift 



space distortions (Jennings et al. 2010 MaruUi et al. 
20Tll. 



3. DATA 

For the analysis presented here we have used the fol- 
lowing publicly available data sets: 

CMBwMAp: The CMB as observed by Wilkinson Mi- 
crowave Anisotrop y Probe (WMAP) fr om the 7 
year data releas^ ( [Komatsu et al.|[20lT ). 



CMBspt: The CMB a s observed by the South Pole 
Telescope^ (SPT) ( [Keisler et al.||201l[ ). 



Hq: a Gaussian prior of Hq = 73.8 ± 2.4 kms"^ Mpc ^ 
on th e Hub ble parameter value today from |Riess| 



et ah (20111. 



H{z): The expansion history of the Universe as measured 
with passively evolving galaxies at z < 1.75. The 



2' 



camb . info 



Available trom |htfr >: //lambda. gsfc ■ nasa.gov 



Available 

[kelslerll/ 



from 



http : //pole .uchlcago ■ edu/publlc/data/| 



combined sample of Simon et a l. (2005); Stern et al 
J^OIO); Moresco et al.| ( |2012b[ ) as given in [Morescq 
et al.i (|2012a|). The expansion history is consistent 



et al.| ( |2012a| . The expansion history 
with H [z) as measured from t he Alcock-Paczy nski 



test with WiggleZ galaxies ( [Blake et al.||2012 |. 



BAO: The BAG scale measured at z = 0.106 {rs/Dy = 
0.336 ± 0.015) fr om the Six Degree F ield Galaxy 
Survey (6dFGS) ( [Beutler et aL||2011[ ), the recon- 
structed value at z = 0.35 [rs/L>v = 0.1126 ± 
0.0022) from Sloan Digital Sky Survey (SPSS) 

Luminous R ed Galaxies (LRG) (Padmanabhan 

et aL||2012l), and z = 0.57 (r,/£)V = 0.0732 ± 
iXD012) from Baryon Oscillation Spectroscopic Sur- 
vey (BOSS) (Anderson et al.,2012j . The measure- 
ments of TsfUv are independent of the assumed 
fiducial values of A^cff (as they cancel out). 

SN: Type la supernovae as o bserved by the SuperNova 



et al.||2011[|Conley et al 



Legacy Survey (SNLS) (Guy et al. 2010 Sullivan 



P(k): The WiggleZ Dark Ener gy Survey power spec- 
trunj^as described in Parkinson et al. (2012). 



3.1. Comments on data 



ibinations 



WMAP and SPT are highly complementary observa- 
tions of the CMB. WMAP is a full sky survey includ- 
ing the largest scales (multipole moments 1 = 2 — 1200) 
while SPT is a high resolution small scale survey {I = 
650 — 2999) very sensitive to effects of changing iVoff (see 
the right panel of Fig. [T]). 

The large scale structure can be included either in the 
form of the BAO scale or the full power spectrum shape. 
To be conservative we do not include both the BAO dis- 
tance measurement and the power spectrum shape data 
from the same survey since we do not know the corre- 
lation between the two. For our main result we use the 
BAO scale from 6dFGS, SDSS and BOSS, and the full 
WiggleZ power spectrum. We have tested that the re- 
sults from this setup are very similar to using the full 
power spectrum shape from both SDSS and WiggleZ. 
This is because the constraints are limited by the deter- 
mination of the ACDM parameters rather than by the 
observational determination of the slope of the galaxy 
power spectrum. The precise determination of the BAO 
scale strengthens the constraints on the ACDM param- 
eters, which improves the ^ rrii, and A'cff constraints by 
breaking degeneracies. Therefore including some BAO 
data is more important than including the full power 
spectrum shape for all LSS data. 

A recent analysis based re-calibratio n of the cepheids 
found Hq = 74.3 ± 2.1 kms~i Mpc~^ ( [Freedman et al. 
2012 1 which is slightly larger than the value adopted here 
driving the value of A^cff upwards. However, the tight- 
est constraints presented here does not rely on the Hq 
measurement but the BAO and SN to constrain the ex- 
pansion history. 

For the SN we use the SNLS data which is a consis- 
tent sample of SN observed with a single telescope and 
conseque ntly free from some o f the ligthcurve fitting is- 
sues that Joudaki et al. (2012) pointed out between the 
Union2 and SDSS SN dafasets. 



Available from http://smp.uq.edu.au/wigglez-data 
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Figure 2. The 68% and 95% CL contours in the N^ff-'^rrii, pa- 
rameter space of fitting a ACDM-I- mi,+NaS model to a selection 
of data combinations. 

4. METHOD 

We fit the data to a standard flat ACDM cosmol- 
ogy with the following parameters: the physical baryon 
density (Vly^h?)^ the physical cold dark matter density 
(r^dm/i^) the Hubble parameter at z = (i?o) the optical 
depth to reionisation (r), the amplitude of the primordial 
density fluctuations (As), and the primordial power spec- 
trum index {us). In addition we vary the effective num- 
ber of relativistic degrees of freedom [N^s — N^ + ANcs) 
and the sum of neutrino masses rrii, = N^rrii,) where 
Ni, is the number of massive neutrinos usually taken to 
be Ni, = 3.046 and the 0.046 accounts for the neutrino 
energies that arise due to the residual heating provided 
by the e"'"e~ -annihilat ions after the neutrino decoupling 



( [Mangano et al.||2005[ ) 

We sample the parameter space usin g the Markov 
Chain Monte Carlo sampleij^ CosmoMC ( Lewis & Bri 



die 



2002D. Details on the specific modules can be found 
'arkinson et al. ( |2012[ ) and on the WiggleZ website^. 
'Tie probability distributions are calculated by 
marginalising over the density of points in the chains. 
Because the distributions can be non-Gaussian, the 2a 
uncertainties can be larger than the 2x la uncertainties. 

5. RESULTS AND DISCUSSION 

From fitting ACBM+J2m^+NcS to CMBwmap + 
CMBsPT + WiggleZ + H(z) + BAG + SN we obtain 
iVeff = 3.58t°:}^(68%CL)iS]j5(95%CL) and Y.^^ < 
0.60eV(95% CL) shown in Fig. [2] and Fig. [3| This not 
represents a preference for more than three neutrino 
species at 95% confidence. 

5.1. Comparison to ACDM 

Fig. Hlshows the resulting contours from fitting 
ACDM+^m^-t-A^eff to CMBwmap + CMBspt + Wig- 
gleZ + H{z) + BAG -I- SN (sohd black), compared to fit- 
ting a pure ACDM model (dashed red), AGBM+J^m^ 
(dotted blue), and ACDM+A"off (dot-dashed green) to 
the same data. Naturally the constraints are tighter for 
pure ACDM, since there are fewer parameters to con- 
strain, but overall the constraints are consistent for all 

® http://cosmologist.info/ 
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Figure 3. The one dimensional probability distributions for A^cff 
and Yl ra^ of fitting a ACDM-f J] 

model to a selection of 
data combinations (same colours as in Fig. [2]l. The grey shaded 
area indicates the probability distribution when correlations are 
neglected by fixing J^m^=0 (left) and Afoff=3.046 (right). 
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Figure 4. The 68% and 95% CL contours for fitting 
ACDM-(-X; m^-l-Afcff (solid black), ACDM-t-J] (dotted blue), 
ACDM-I-Afcff (dot-dashed green), and ACDM (dashed red) mod- 
els to CMBwmap + CMBspt + WiggleZ 4- ll{z) + BAO 4- SN 
data. The resulting contours arc consistent with each other for all 
parameters. 

parameters. Adding ^ to the fit does not change the 
ACDM parameters or uncertainties significantly whereas 
adding iVeff increases the values and uncertainties of iJg. 
Adding both parameters increases the preferred values of 
iJp and Ug but also the uncertainties, so the values remain 
consistent with the pure ACDM case. The shift can be 
understood as follows: Increasing iJo changes the height 
of the peaks in the CMB, as it corresponds to increasing 
the physical matter density. Changing iVeff recovers the 
details of the CMB peaks (because the original ratio of 
matter to radiation is recovered, restoring the details of 
Silk damping), but with too much power on large scales 
due to ISW. This shifts the primordial power spectrum 
from very red (ris — 0.96) to slightly less red (0.98— 1.0). 

5.2. Combinations and degeneracies 

It is clear from the 'Yrn^-Nf.g contours in Fig. [2] that 
unless all data sets are included, the two parameters are 
correlated. The tightest constraints come from the com- 
bination of CMBwmap + CMBspt + WiggleZ + H{z) 

BAG -|- SN. Fig. [5] demonstrates the degeneracies be- 
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Figure 5. Results of fitting a ACDM+ m^+A^cff model to a 
selection of data combinations (same colours as in Fig.[2]|. 

tween J2 & ^cff , and i/o & ^^m- When all the data is 
combined, there is a significant correlation between Hq 
and iVeff, and a mild correlation between ^ and fi„i- 
Lowering Hq shifts iVog back towards the expected value 
of three. 

At the present, large scale structure (e.g. WiggleZ 
and SDSS) does not constrain Ncs uniquely because the 
turnover of the power spectrum is not well measured, 
which leaves some correlation between fim, X]'^:^) and 
N^ff (see FigJTl), but inclusion of the power spectrum 
improves on YZ^f from the small scale suppression. 

The H(z) data does not add much to the neutrino 
mass constraints, but because it constrains the expansion 
history, it improves on the limits on N^s as discussed in 
Moresco et al.l (|2012a|. 



5.3. Nucleosynthesis 

There have been a num ber of recent attempts to derive 
A^eff from BBN alon e (e.g. [Mangano fc Serpico|2011[ [PeF 
tini fc Cooke||20l"2 ) which ail seem to be consistent with 
Aeif = 3. Analys es combining BBN and CMB seem to 



prefe r N >3 (e.g. [Nollett fc Holder [201 l[ |Hamann et al. 
20111. It is concerning that the preference tor l\eS> is 
present in all analyses including CMB, but not preferred 
by BBN alone. This could indicate a systematic error 
in one of the data sets. However, BBN alone relies on 
very few and notoriously difficult measurements of the 
deuterium and helium fractions. 

The difference between the results could also be in- 
terpreted as potential tension between A'off measured at 
two different epochs (BBN and recombination), where we 
naively would expect A'eff to be the same. A temporal 
variation can be explai ned theoretically by a decaying 
particle ( Boehm et al.lj2012 ) but currently there is no 
experimental evidence for the existence of such particle. 

5.4. Beyond ACDM+J2 +Ncf[ 

Other physical effects such as curvature, varying equa- 
tion of state, running of the spectral index etc. could 
mimic Npff> 3 i f not p roperly accounted for in the mod- 



elling. Joudaki (2012) demonstrated that the deviation 



CMBsPT, BAO scale from SDSS and 2dFGS, Hq from 
HST, and SN from Union2. However, for all parameter 
combinations the preferred value of A'off was still above 
three, and only when more than one extra extension of 
the ACDM cosmology (e.g. curvature and varying equa- 
tion of state in addition to N^g and ^m^) were consid- 
ered, did the preferred value of A"cff become consistent 
with three within one standard deviation. We take this 
as an indication that one extension alone does not ex- 
plain the preference for A'off 3. 

5.5. Future 

Measuring the position of the peak of the matter power 
spect rum (the turnove r) would give another handle on 



Noff. Poole et al. 



from Aeff = 3.046 is diminished if allowing for curva- 
ture, varying equation of state, running of the spectral in- 
dex, and/or the helium fraction by fitting to CMBwmap, 



S2012 ) predicts that a Euclid-like galaxy 
e to constrain N^g to approx. 20% in- 
dependently of the CMB. However, the position of the 
turno ver can be degenerat e with neutrino hierarchy ef- 
fects (Wagner et al. 2012). The degeneracies between 
^ TO^ , Aofi, and the hierarchy allows for a measurement 
of either of them if the remaining parameters can be 
"fixed" by independent data. 

If A^cff can be measured with sufficiently high precision, 
eventually it will be possible to measure the thermal dis- 
tortion of the neutrino spectrum (the 0.046). 

Any ncutrino-like behaving particle, including sterile 
neutrinos and axions etc., which decouple early when 
relativistic and become non-relativistic, can mimic the 
effect of the neutrinos. If the value of N^s > 3 is due 
to a particle, the result points towards physics beyond 
the Standard Model, and its existence it will have to be 
confirmed by a particle physics laboratory experiment. 
If no sterile neutrinos are found in laboratory experi- 
ments, the cosmological preference for additional species 
may indicate a lack of understanding of early Universe 
physics. 

6. CONCLUSION 

Due to imperfect measurements of the ACDM parame- 
ters ^ and A'eff are not entirely independent param- 
eters and should be fitted simultaneously in cosmologi- 
cal analyses. Performing such fit of ACDM-|-^ mi^+N^g 
to a combination of cosmological data sets, leads to a 
2cr preference for Aeff> 3.58t^ j^(95%(CL) and ^ < 
0.60 eV, which are currently the strongest constraints on 
A'off from cosmology simultaneously fitting for N^s and 
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